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Abstract 

1^ \ The partition function of the Baxter- Wu model is exactly related 

0^ ' to the generating function of a site-colouring problem on a hexagonal 

'^ ■ lattice. We extend the original Bethe ansatz solution of these models 

^s^ i in order to obtain the eigenspectra of their transfer matrices in finite 

^D '_ geometries and general toroidal boundary conditions. The operator 

content of these models are studied by solving numerically the Bethe- 
ansatz equations and by exploring conformal invariance. Since the 
Cd I eigenspectra are calculated for large lattices, the corrections to finite- 

size scaling are also calculated. 

c . 

o 1 Introduction 

o 
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The Baxter- Wu model is defined on a triangular lattice by the Hamiltonian 

H = -J Yl aiffjak, (1) 

<ijk> 

where the sum extends over the elementary triangles and ctj = ±1 are Ising 
variables located at the sites. This model is self-dual ill with the same 



critical temperature as that of the Ising model on a square lattice, and was 
solved exactly in its thermodynamic limit by Baxter and Wu . Its leading 
exponents a = 2/3, /i = 2/3 and rj = 1/4, are the same as those of the 
4-state Potts model 0, §, ^. Due to this and the fact that both models 
have a fourfold degenerate ground state, it was conjectured that they share 
the same universality class of critical behaviour. However from numerical 
studies of these models on a finite lattice it is well known that both models 
show different corrections to finite-size scaling. While in the Potts models 
[i? 07 i) m these corrections are governed by a marginal operator, producing 
logarithmic corrections with the system size, this is not the case in the Baxter- 



Wu model [|T^, |T^. This raises the question of knowing which operator 
governs these corrections in the Baxter- Wu model. 

With the developments of conformal invariance applied to critical phe- 
nomena [|r^, two models are considered in the same universality class of 
critical behaviour only if they have the same operator content, not only the 
leading critical exponents. The operator content of the 4-state Potts model 
was already conjectured from finite-size studies in its Hamiltonian formula- 



tion |T^, and can be obtained by a Z{2) orbifold of the Gaussian model (see 



n] for a review). 

In this letter, by exploiting the conformal invariance at the critical point, 
we report on our calculation of the operator content of the Baxter- Wu model. 
In order to do this calculation we have to generalize the original Bethe ansatz 
solution of the model, since this solution does not give the complete eigen- 
spectrum of the associated transfer matrix T. The conformal anomaly c and 
anomalous dimensions (xi,X2, ...) are obtained in a standard way from the 
finite-size behaviour of the eigenspectra of the associated transfer matrix, at 
the critical temperature. If we write T = exp(— if ), then in a strip of width 
L with periodic boundary conditions the ground-state energy, Eo{L), of H 



behaves for large L as ||T5| 



M^ = ,_|^+„(i-.), (2) 

where eoo is the ground-state energy, per site, in the bulk limit. Moreover, 
for each operator O^ with dimension Xa there exists a tower of states in the 
spectrum of H with eigenenergies given by |T^ 



E^AL) =Eo + —^{x^ + m + m') + o{L~^), (3) 



where m,m' = 0,1,2,... . The factor Vg appearing in (^ and (^ is the 
sound velocity and has unit value for isotropic square lattices. The higher 
eigenvalues of T can be calculated directly by numerical diagonalization. 
However since T is not sparse and has dimension 2^, for a horizontal width 
L, we cannot compute its eigenspectra by numerical diagonalization methods 
for L > Lq r^ 26. Instead of a direct calculation we relate this problem 
to a site- colouring problem (SCP) on a hexagonal lattice, which can be 
solved by the Bethe ansatz. Following P the partition function Z^^^, of the 
Baxter- Wu model on a periodic triangular lattice with L (N) rows (columns) 
in the horizontal (vertical) direction is related to the partition function, or 
generating function, Z^fJ'j^ of a SCP on a hexagonal lattice with M = 2L/3 
(N) rows (columns) in the horizontal (vertical) direction. In the limit A^ -^ 
oo, 

'vBW _ ySCP {A\ 

^LxN — ^My.N- V^) 

In figure 1 the Baxter- Wu model is defined on the triangular lattice formed 
by connecting the points o, o and A. The related SCP is defined on the 
hexagonal lattice formed by the points A and O, connected by continuous 
lines. The configurations in the SCP are defined by attaching at the points of 
the hexagonal lattice (i, j) (i = 1, . . . , M; j = 1, . . . , A^) site-colour variables 
Cjj- = 1, 2, . . . , 8, satisfying the constraint that any two nearest-neighbour 
colours must differ by 1 or —1. The partition function Z^'-^^ given in (^) 
is obtained by adding all the colour configurations with weights given by 
the product of the fugacities Zj{j = 1,2, ...8) of each colour in the lattice 
configuration. These fugacities are given by p[ 

zi = zs = Z5 = zr = 2 sinh(4/3 J) ,^. 

Z2^ = Z4 = Zq^ = Zg = sinh(2/?J) = t. ^ ' 

The critical point of the Baxter- Wu model and of the SCP is given by the 
self-dual point t = t^ = I. If we write T = exp{—H) for both models and 
since Z = Tr{T^), the relation (^) implies 

Tr(e-^^"")=Tr(e-^^M''") (6) 

Its is important to observe that although H^^ and H^f^ have the same 
dimension 2^ they may have different eigenvalues. 

In the SCP we say we have a dislocation in a given link wherever the 
colour in its right end is smaller than the left one. The number n of disloca- 
tions in a given row is a conserved quantity. Consequently the Hilbert space 



associated to T^'~^^ or H^^^ is separated into block disjoint sectors labelled 
by the values of n. For a periodic hexagonal lattice of width M (even) the pos- 
sible values of n are even and given by M, M±4, M±8, . . . , M±4 int{M/4). 
Due to (§) the SCP has also an additional Z{2) symmetry (eigenvalues 
e = ±1), since adding 4 (modulo 8) to all colours in a given configuration 
does not change its weight in the partition function. The Bethe ansatz solu- 
tion presented by Baxter and Wu [Q only gives part of the eigenspectrum of 
ffSCP^ since they considered only eigenstates which are even under this sym- 
metry. We generalized ^^ their solution in order to obtain the missing odd 



eigenvectors. For the sake of brevity we only present here the Bethe ansatz 
equations. The eigenvalues En of H^^^ in the sector with n dislocations 
are given by 

/\/f " 

Et^ = -^ ln(16t^(l + t^)) - E(e^^ - ^kf^\ (7) 



where 



ej'^'^ = 1/2 In (xj + Sj^xj - l) ; Xj = cos{2kf'^) + t + 1/t, (8) 



with 1 = si = S2 = . . . = Sn-i = -Sn^i+1 = . . . = -s„, and / = 0, 1, ...,n 
fixed. The quasi- momenta {kj } are obtained by solving the equations 

/ (s)\ ^ " /cosh(e^'^ + iA;("''))\ 

exp .Mfcf-) =-Vin , M^-X.) h ^■ = l'2,...,n, (9) 

p=i \cosh(ep ' +iky')J 

where e = ±1. The value e = 1 gives the part of the eigenspectrum which is 
even under the Z{2) symmetry and was derived in 0, while e = — 1 gives the 
odd part of the eigenspectrum. Strictly speaking this is a conjecture since 
the completeness of the Bethe ansatz solutions is always a difficult question. 

We have studied extensively these equations numerically at the critical 
point t = tc = 1, for general values of n, e and / and for lattice sizes up to 
M ~ 200. For example, in the eigensector where n = M, by setting / = 
we obtain the energies of the ground state and first excited state by choosing 
e = 1 and e = — 1, respectively. 

Let us calculate the conformal anomaly by using @. The bulk energy 
e^^ = — I In 6 can be obtained from the solution in the bulk limit and the 



sound velocity Vs = "\/3/3, can be inferred from (^) and an overall analysis 
of the dimensions appearing in the model. Using these values in (0) we 
obtain the finite-size sequence estimators c{M) for the conformal anomaly, 
presented in table 1. As expected the conformal anomaly is c = 1, like the 
4-state Potts model. The relation (|^) does not imply that in a finite lattice 
the eigenergies of H^^ and H^f^ are the same. In fact this is the case. 
Fortunately by comparing their eigenspectra by direct calculations on small 
lattices we verify that many of the lower energies, including the ground-state 
energy, are exactly the same. Consequently by using the bulk limit value 
e^ = — |ln6 we obtain the same sequence shown in table 1. The sound 
velocity which comes from our analysis is now Vg = a/3/2 and the conformal 
anomaly has the expected value c = 1. 

The large-L behaviour of the energies of excited states will give the op- 
erator content of the models. Using (^, the finite-size sequences obtained 
for some dimensions associated to H^^^ are shown in table 2. The estima- 
tors x'j{M — n,l) are the j*^ lowest eigenenergy obtained by solving (0)-(§) 
with the values n, e and I. As a result of an extensive calculation of the 
eigenspectra of H^'-^^, obtained by direct diagonalization on small lattices 
and by solving (0)-® for large lattices, we arrive at the following conjecture. 
Namely the dimensions of primary operators in a given sector labelled by 
n = M + 4p are given by 
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a;p,, = ^(V + ^), g = 0,±l,±2,... , (10) 

where p = 0, ±1, ±2, . . . for the periodic lattice. The number of descendants, 
with dimensions Xp^q + m + m' (m, m' G Z) is given by the product of two 
independent Kac-Moody characters. The dimensions ( p!0D are similar, in a 
Gaussian model |18[ with compactification radius equal 2, to the dimensions 
of an operator with vorticity p and spin- wave number q. The Gaussian model 
at this radius corresponds to the continuum limit Kosterlitz-Thouless point 
of the X-Y model in the torus [[19], ^M. This implies that the SCP belongs 



to the same universality class as the Kosterlitz-Thouless phase transition of 
the X-Y model. 

We have also studied the SCP with more general toroidal boundary con- 
ditions which preserves the same symmetries as the periodic case ( n and e 
are good quantum numbers). Those boundary conditions are obtained by 
imposing to each row (j = 1, 2, . . .) of a colour configuration in the SCP the 



constraint cm+i,j = Cij + k, where k = 0, 2, 4 or 6. The periodic case is 
obtained when k = 0. The possible values of n are now given hj n = M + 4p 
where p = j — \ (j = 0, ±1, ±2, . . . ). In this case k = 2 and A; = 6 we were 
not able to apply the Bethe ansatz for arbitrary temperatures, but only at 
the critical temperature t = tc = 1 [0. The Bethe ansatz equations turn 



out to be the same as in (0)-® but now the values of n depend on the 
boundary condition. If in @ we take Eq{M) as the ground-state energy of 
the periodic lattice {k = 0) our numerical solutions of Eqs. (0)-® also give 
the dimensions (|10D but with p = j — \ (j = 0, ±1, ±2, . . .). 

Let us return to the Baxter- Wu model. In this case by comparing the 
eigenspectra of H^^ and H''^^^, obtained by a direct diagonalization on 
small lattices, we verify that many of the dimensions Xp^q appearing in (^^) 
are absent. For example the energies producing the estimators in the 2"*^ 
and 8th columns of table 2 only appear in H^^^ . Following for large lattices 
the energies which are exactly related in both models we verified that the 
lower dimensions in the Baxter- Wu model are given by x = 0,|,|,|,..., 
and appear with degeneracy dx = 1,3,1,9,..., respectively. These results 
supplemented with the global eigenspectrum calculated for small systems, 
indicate that the operator content of the Baxter- Wu model is the same as 



that of the 4-state Potts model [0 and is given in terms of a Z{2) orbifold 
[|T^ of the Gaussian model. It is interesting to remark that while in the 
SCP the operator content is given in terms of characters of the Kac-Moody 
algebra in the Baxter- Wu model the characters are those of the Virasoro 
algebra. Since the exact integrable SCP and the Baxter- Wu model belong 
to different universality classes some care has to be taken when we import 
exact results from the SCP to the Baxter- Wu model. 

We also studied the Baxter- Wu model with more general toroidal bound- 
ary conditions which preserve its Z{2) ® Z{2) symmetry. We observe numer- 
ically that the eigenenergies which appear in this case can also be obtained 
from the eigenspectrum of the SCP with the toroidal boundary conditions 
we considered {k = 0,2,4,6). Calculating the corresponding dimensions for 



large lattices we obtain the same dimensions reported in |TB[ for the 4-state 
Potts model. These results imply that both models are indeed in the same 
universality class, being governed at the critical point and arbitrary toroidal 
geometry by the same conformal field theory. 

Since we calculate eigenenergies of H^^ and H^^^ for large lattices we 
can now also calculate the corrections to finite-size scaling for both models. 
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Consider the lowest eigenenergy Ea, associated to an operator with dimension 
Xa- From (|]) the correction Ra{M) associated to this level is given by 

E^ = Meoo + ^ (^" - Y^ + Rc.{M)j . (11) 



According to conformal invariance |12, H i?„ should behave as 



floW-Eirs^rj + ETfai:;^. d^) 



' 7:7 

where {x^} are the non-relevant dimensions (x^ > 2) associated to the op- 
erators governing the finite-size corrections. In the 4-state Potts model the 
lowest dimension x^ in this set is associated to a marginal operator (x^ = 2), 
and the corrections have a logarithmic behaviour with the system size. In 
table 3 we show our estimators for the dimension x^ of the dominant correc- 
tion for some eigenergies Ea, with corresponding dimension Xa, in the SCP 
and the Baxter- Wu model. In all these cases we clearly see that x^ = 4 
indicating that the corrections are integer powers. Rather than the 4-state 
Potts model these corrections are like those of the Ising model. This explains 
why the finite-size studies of the Baxter- Wu model have good convergence, 
in contrast with the 4-state Potts model. 

These results show that although the 4-state Potts model and the Baxter- 
Wu model share the same universality class of critical behaviour, having the 
same operator content, the finite-size effects correspond to different pertur- 
bations of the fixed point of the renormalization group. The SCP belongs 
to another universality class. This implies that not all exact results derived 
from the SCP should be translated to the Baxter- Wu model. 
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Figure and Table Captions 

Figure 1 - The Baxter- Wu model is defined on the triangular lattice formed 
by the points o, o and A. The site-colouring problem (SCP) is defined on 
the hexagonal lattice formed by the points A and O, connected by continuous 
lines. 



Table 1 - Conformal anomaly estimators cm, as a function of M, for the SCP 
and Baxter- Wu model. 



Table 2 - Scaling dimensions estimators x^(M — n,l), as a function of the 
lattice size M, for some eigenenergies. These energies are the j^^ lowest 
energy obtained by solving (0)-(|D with values n, e and /. 



Table 3 - Estimators of the exponent x^ in (12) of the dominant correction 
of some eigenenergies E^ with dimension Xa- 
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Table 1 



M 


Cm 


6 


0.996590995 


10 


0.998910268 


50 


0.999959561 


100 


0.999989915 


150 


0.999995519 


200 


0.999997480 



Table 2 



M 


a:r(0,0) 


X2-(l,0) 


a;+(0,l) 


a;2+(2,0) 


x+(0,l) 


6 


0.12502803 


0.24896741 


0.50626226 


0.62613504 


0.98648357 


10 


0.12501702 


0.24959771 


0.50215317 


0.62548322 


0.99698967 


50 


0.12500083 


0.24998323 


0.50008406 


0.62502093 


0.99992101 


100 


0.12500021 


0.24999580 


0.50002099 


0.62500524 


0.99998057 


150 


0.12500009 


0.24999813 


0.50000933 


0.62500233 


0.99999139 


200 


0.12500005 


0.24999895 


0.50000525 


0.62500131 


0.99999516 
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Table 3 



M 


Xa = 


Xa = 0.125 


Xa = 0.25 


Xa = 0.5 


Xa = 0.625 


20 


4.0836238 


4.7672008 


3.953517 


4.0274781 


3.9063776 


50 


4.0466061 


4.8743240 


3.990329 


4.0576819 


3.9816761 


100 


4.0035431 


4.9917793 


3.998189 


4.0010866 


3.9965769 


150 


4.0012060 


4.9977042 


3.999406 


4.0003372 


3.9989165 


200 


4.0005332 


4.9991458 


3.999680 


4.0001660 


3.9994630 



12 



